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Abstract 

We report on the study of the mass-radius (M-R) relation and the 
radial oscillations of proto strange stars. For the quark matter we have 
employed the well known density dependent quark mass model and 
its very recent modification, the temperature and density dependent 
quark mass model. We find that the maximum mass the star can 
support increases significantly with the temperature of the star in this 
model which implies that transition to a black hole at the early stage 
of formation of the star is inhibited. As for the neutrinos, we find, 
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contrary to the expectation that the M-R and oscillation frequencies 
are almost independent of the neutrino chemical potentials. 
Subject headings: Strange stars - Oscillations 
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1 Introduction 



It is believed at present that at high densities (p > 3po) the 
hadronic matter undergoes a phase transition to the unconfined 
state of quarks and gluons which is composed of roughly equal 
number of u, d and s quarks together with a gas of electrons and 
muons necessary to maintain charge neutrality. This strange 
quark matter (SQM) may be the true ground state of matter 
(Witten 1984; Farhi & Jaffe 1984; see also e.g. Madsen 1999, 
for a recent review of the physics and astrophysics of strange 
quark matter and strange stars). Since then, various implica- 
tions of SQM have been explored in astrophysics and high energy 
physics. Glendenning (1997) has conjectured the nucleation of 
SQM bubbles and most of the currently believed neutron stars 
may be strange stars. 

Pulsars can be modeled as neutron stars or strange stars. It is 
very difficult to distinguish observationally between the two of 
them. Recently three ways have been proposed for identifying 
a strange star: 

1. A hot strange star or a proto strange star (PSS) may rotate 
more rapidly since its higher bulk viscosity (Wang and Lu 1984; 
Madsen 1992; Anand et al. 2000a) can effectively damp away 
r-mode instability (Madsen 1998). 
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2. Mass-Radius (M-R) relations of strange stars (M oc R 3 ) is 
quite different from that of neutron stars (Li et al. 1999). 

3. It is possible to distinguish "bare" polar strange stars from 
neutron stars via pulsar magnetospheric and polar radiation (Xu 
et al. 1999; Xu et al. 2001). A newly born star, be it a neutron 
star or strange star passes through various stages of evolution. 
The duration of these stages is essentially determined by neu- 
trino diffusion time scales (Prakash et al. 1997). Roughly the 
diffusion time scale is proportional to R 2 X where R is the star 
radius and A the neutrino mean free path. Both the M-R re- 
lation and mean free path depend upon the equation of state 
(EOS) and the composition of the star matter. The behaviour 
of the maximum mass as a function of the temperature and neu- 
trino trapping is also of practical importance to the process of 
conversion of the star into a black hole and the stage at which 
this transition occurs, if at all. 

The formation of proto neutron stars (PNS) is relatively well un- 
derstood (for a recent review see e.g. Prakash 1997; Glendenning 
1997). But due to the absence of a full theory determining the 
conditions at which quark matter phase transition occurs and 
a lack of detailed understanding of the complex burning pro- 
cess of neutron matter into strange matter, very little is known 
about how a PNS converts to PSS during the supernova explo- 
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sion. However some studies of the transition of nuclear matter 
to SQM taking into account the two step phase transition viz 
neutron matter — > ud — > uds have been carried out by numer- 
ous authors (Dai et al. 1995; Anand et al. 1997; Lugones and 
Benvenuto 1998,1999). 

In the first step the neutron matter deconfinement occurs on 
a stronge interactions time scale ~ 10~ 23 s; after that chemical 
equilibrium of the deconfmed matter takes place on a weak in- 
teraction time scale ~ 10~ 8 s. Additional neutrinos and energy 
are produced in the second step. Recently Gupta et al. (2002) 
have studied the phase transition in a PNS from nuclear matter 
to two flavour quark matter. 

Most of the calculations for the phase transition and the prop- 
erties of strange stars have been done in the MIT Bag model. 
There is an alternative and equally attractive approach based 
on the density dependent quark mass (DDQM) model originally 
proposed by Fowler, Raha & Weiner (1981). Many aspects of 
SQM like radial oscillations of both rotating and non rotating 
quark stars, bulk viscosity etc. have been studied by various au- 
thors in this model (Anand et al. 2000a; Anand at al. 2000b). 
At zero and very low temperatures the model works reasonably 
well and gives results for various quantities like the M-R relation 
and viscosity etc. which are qualitatively similar to those ob- 
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tained with bag model. Very recently Zhang & Su (2002) have 
pointed out that the DDQM model suffers from a basic draw- 
back that it can not be used to explain the process of quark 
deconfinement phase transition as the quark confinement is per- 
manent in this model. They have suggested a temperature and 
density dependent quark mass (TDDQM) model to describe the 
bulk SQM. In this model the mass of quarks not only depends 
on baryon density (as in DDQM model) but also on the tem- 
perature. In DDQM model the masses of u, d and s quarks are 
taken to be (Fowler, Raha & Weiner 1981; Lugones and Ben- 
venuto 1998) 

C 

m u = m d = - — , 
3n B 

C 

m s = m so + - — (1) 
3n B 

In TDDQM model guided by Friedman-Lee soliton bag model 
(Lee 1981), C depends on temperature and can be parameterized 
as 

C(T) = C [l-(^ 2 ]; < T < T c 

-L c 

= 0; T > T c (2) 

where T c is the critical temperature of deconfinement. At T = 
0, it of course reduces to the usual DDQM model. One system 
where the results can be expected to be affected significantly 



by the inclusion of this temperature dependence of the quark 
masses is the PSS which are born with a rather high temperature 
of ~ 60 MeV or more and large lepton excess. Motivated by 
these factors we have studied in this paper the M-R relations and 
radial oscillations of PSS during the first few (~ 20) seconds of 
its existence when it has a high temperature and large number 
of electron (and possibly muon) neutrinos. In section 2, we 
present the formalism and the thermodynamics used. Section 
3 deals with M-R relation and radial oscillations of PSS and 
section 4 is devoted to results and discussion. 

2 Formalism and Thermodynamics 

We consider PSS as a system of uds quarks along with electron, 
muon and their neutrinos. The thermodynamical potential for 
the system is given by Q = E where 

Qi = -0/^[log(l+exp(-^ -m])) 

+ log(l + exp(-/p7i + &]))] (3) 

and 

Ei = \jp 2 + m? (4) 

Qi= 6 (3 colour x 2 spin) for u, d and s quark, two for e and \i and 
one for n„ e and {i v . For quarks, the expressions for the masses in 
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the TDDQM model are given by equation (1) with C replaced by 
C(T) of equation (2). The thermodynamic quantities of interest 
like the energy density, p, the pressure, P and various number 
densities rii can be obtained from the thermodynamical potential 

d{n t /n B ) dtt t dQidrrii 

P * = — ^TTi — r = n B^ n i = n B^ — o (5) 

0(1 /riB) oub orjiidnB 

The additional term ns§^ arises because of the dependence of 

the quark masses on the baryon density. 

pi = n { + piUi - T-^ (6) 

and number density rii by 

ni = —z— (7) 

Using eqn (3) for the thermodynamical potential we get, for the 
quarks 

p = \ 9i_ f P 4 dp[fj(T) + fj(T)} 
32tt 2 J E % 
C g l m l r p 2 dp[fj(T) + fj(T)} 
3n B 2tt 2 J E % K ' 

9i jp 2 Eidp(fi(T) + UT)) 



Pi = 



2tt 2 



1 9 l m l C T 2 f p 2 dp[fj(T) + ]j(T)] 

3 7r 2 n B T 2 J E % 1 ] 



9i 



I p 2 dp{f,(T) - UT)\ (10) 



Here fi(T) and fi(T) are the fermi-dirac distribution functions 
for particles and antiparticles respectively: 

MT) = (exp([ v / ^ + m2 - ^]/T) + l)" 1 (11) 

UT) = (exp([v^ + m? + ^]/T) + l)" 1 (12) 

The pressure and energy of electrons, muons and their neutrinos 
is given by the same expressions (8) and (9), except that the 
second term does not contribute. The total Pressure, P and 
energy density, p are the sum of contributions from each species 
of particles: 

P = EP, + -Lr 2 T 4 (13) 
45 

P = Epi + ^T 4 (14) 

where the last term is the contribution of the gluons. 
Having derived in detail the thermodynamics with variable par- 
ticle masses, we now apply it to the investigation of PSS. Since 
PSS is a mixture of u, d, s quarks, electrons, muons and their 
neutrinos with (3 equilibrium maintained by the weak interac- 
tions, 

d,s <-> w + e + z/ e , s + w <-> u + d 
u + fJL + 

the various chemical potentials \i{ (i =u, d, s, e, v e and v^) satisfy 

fid = Ps, Pd = Pu + He — Hv e 
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H - AS = Ve- /S (15) 

Charge neutrality of the system demands 
2 11 

^n u - -rid ~ ^n s -n e -n^ = (16) 

The total baryon density is given by 

n B = 7^(n u + n d + n s ) (17) 

Equations (15 - 17) can be solved self consistently for any given 
baryon density, temperature and the two neutrino chemical po- 
tentials. Once these equations have been solved the pressure, P, 
and energy density, p, can be obtained from equation (13) and 
(14) and these give the required EOS, P(p). 



3 Radial Pulsations of PSS 

The radial pulsations of a non rotating star in general relativity 
were first studied by Chandrasekhar (1964). The metric used is 

ds 1 = -e v c 2 dt 2 + e x dr 2 + r 2 {d6 2 + sin 2 6d(j) 2 ) (18) 

In hydrostatic equilibrium, the structure of the star is described 
by the Tolman-Openheimer-Volkoff equations 

dm 9 

— = 4?rr 2 p 19 
dr 
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dP ^ -G(P + pc 2 ) (m + 47rr 3 P(T 2 ) 

~dr~ ~ c 2 r 2 (l - 2§M) ( } 

di/ _ 2GM(1 + 47rr 3 Pmc- 2 ) 

^ " CV(1 - ( } 

For given equation of state P(p), equations (19-21) can be nu- 
merically integrated outwards starting from r = with a given 
central density p c and the corresponding pressure P c (p c )- The 
equations are solved upto the point where P = which defines 
the radius, R, of the star. The value of m(R) = M is the gravi- 
tational mass of the star. 

We next look at the radial oscillations of the PNS. If 5r is the 
radial displacement, 

« = -, C = r 2 e-/ 2 J (22) 
r 

and projecting the time dependence as exp(iat), one gets the 
equation governing radial adiabatic oscillations (Chandrasekhar 
1964; Datta et al 1998 and Anand et al 2000a) 



where 



e i/-A 



F = —^—-2^ P ) (24) 
P + pc 1 



G = 



P + pc 

H = 



~ x \ P(x+il/)+ d Jm_i (TP) 



2 dr r 



e v - x [4dP 87rGe x P(P + pc 2 ) 1 ,dP s2 



P + pc 2 [r dr c 4 P + pc 2 dr 
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(25) 
(26) 



A is related to the metric function through 



= (27) 



rc 2 



rP = -^-TZ ( 28 ) 



and 

P + pc 2 dP 
c 2 dp 

The boundary conditions to solve the pulsations equations (23)- 
(28) are 

((r = 0) = (29) 

5P(r = R)=0 (30) 

Equation (23) subject to boundary conditions (29) and (30) rep- 
resent a Sturm-Liouville eigenvalue problem for a 2 . The radial 
oscillations equations (23) to (26) are totally model indepen- 
dent and are infact the same whether one considers a neutron 
star, quark star, PNS or PSS. The nature of radial oscillations 
depends only on the equation of state used. 



4 Results and discussion 

In the TDDQM model there are three free parameters C , m so 
and T c . C Q and m so are fixed by the requirement that at T = 
uds should be absolutely stable, i.e. Energy/baryon < 930 
MeV while for ud, E/n B > 940 MeV. This gives us a stability 
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window for C Q and m so . We have chosen two sets, (i) C = 185 
MeV/m- 3 and m so = 150 MeV, and (ii) C Q = 210 MeV/m" 3 
and m so = 100 MeV, which both lie within the stability win- 
dow. As explained in the introduction, following Zhang and Su 
(2002) we have chosen T c = 170 MeV. The chemical potentials 
are evaluated self consistently for a given baryon density, tem- 
perature and neutrino chemical potentials by demanding charge 
neutrality and j3 equilibrium conditions for a given set of values 
of C and m so . This provide us a profile of pressure P and energy 
density p. Using the EOS i.e. p, P profile, we then solve the 
stellar structure equations that provide PSS configurations. The 
numerical integration was done by using Runge-Kutta integra- 
tion procedure to obtain p,m and v as functions of r. We then 
calculate (rP) directly from the EOS for all densities by using 
quadratic difference formula for derivative (^r). The equations 
(23)-(28) were also solved by Runge-Kutta integration method. 
By changing the number of input points, it was estimated that 
a 2 is accurate to one part in 10 3 . 

In fig. 1, we plot the gravitational mass M in terms of solar 
mass vs radius R in Km for p Ve = p v = MeV. Curves a, b, 
c, and d correspond to the temperatures T = 0, 20, 40 and 60 
MeV respectively in fig. 1 to fig. 4. It is seen that the maxi- 
mum mass and the radius increases with the temperature. The 
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curve a) for T = corresponds to the normal strange star which 
has been studied in the same model by Benvenuto and Lugones 
(1998) and Anand et al. (2000a). The present results are differ- 
ent from those obtained by the above authors, since, as has been 
pointed out by Peng et al. (2000) and Zhang k Su (2002), the 
expression for the energy density used by these authors is dif- 
ferent and inconsistent in that the lowest energy state does not 
correspond to zero pressure. Thus the present results must be 
considered more realistic in the DDQM or the TDDQM models. 
Fig. 2 to fig. 4 show the graphs of M vs. R for three other neutrino 
chemical potentials (fi, Ve , ^): (0,200); (200,0) and (200,200) 
MeV respectively. The general features are the same as those in 
fig. 1 . It seems that M-R relations are not too sensitive to the 
neutrino chemical potentials, as against PNS whose structure is 
expected to depend quite significantly on the presence of neu- 
trinos (Prakash et al. 1997). For this reason for the rest of the 
investigations we consider \i Ve = n v = 200 MeV only. 
Fig. 5 deals with the variation of maximum mass with the tem- 
perature. We find that the maximum mass increases signifi- 
cantly with increase in temperature implying that PSS supports 
higher mass when newly born. This, according to Prakash et al. 
(2000) implies that transition to black hole can not take place 
very early; if it takes place at all, it could happen later on ther- 
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mal diffusion time scale ~ 10 s. In fig. 6 we show the effect of 
the strange quark mass on our calculation. We chose two sets 
of C and m so (185, 150); (210, 100) which both correspond to 
the same minimum energy per baryon, viz ~ 924 MeV. We see 
that the maximum mass of the PSS vary by about 10% over this 
range of strange quark masses. In this figure we also show the 
effect of inclusion of the temperature dependence of the quark 
masses. Curve a) corresponds to the TDDQM model whereas 
curve c) corresponds to DDQM model. Inclusion of the temper- 
ature dependence of the quark masses increases the maximum 
mass as well as the radius by about 5-10%. 
In fig. 7, we plot frequency of radial pulsations for fundamental 
(n = 0) and first excited (n = 1) mode vs mass of PSS for ji Ve 
= fi Vii = 200 MeV and temperatures T = and 60 MeV. We 
notice that at any temperature the fundamental mode frequency 
is exactly equal to zero for maximum mass of the star for that 
temperature implying that the numerical code used by us is in 
order. At higher temperatures the frequencies of cold and very 
hot PSS are almost same implying that the hotness of the star 
does not influence much the radial oscillations of the star. In 
fig. 8 we show the effect of strange quark mass and temperature 
dependence on the radial frequencies for n=0 and n=l modes. 
Curves a) and c) correspond to the two sets of C and m so (185, 
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150); (210, 100) respectively for n=0 and (d,f) for n=l. We 
notice that the effect is rather small but relatively more pro- 
nounced for n=l. Same is true for the effect of temperature 
dependence of the quark mass which can be obtained by com- 
parision of curves (a, b) for n=0 and (d, e) for n=l. 

To conclude, we find that the correct thermodynamic treatment 
(Peng et al. 2000) of density dependent masses changes the 
results significantly. The PSS being at a higher temperature 
comparred to normal SS has a much higher maximum mass and 
can thus support a star of higher mass, implying that transition 
to black hole at a very early stage is inhibited. However, con- 
trary to the usual expectation neutrino chemical potentials do 
not seem to have much effect on M-R and the radial pulsations 
of the PNS. Thus the structure depends upon the temperature 
with which the PSS is born and not on the lepton fraction. Since 
the magnetic field has a considerable effect on the possibility of 
phase transition (Gupta et al. 2002), it will be interesting to 
study its effect on the structure of PSS. This work is in progress 
and will be reported elsewhere. 
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Figure captions 



Figure 1. Radius R in km vs. mass M/M© (where M Q is 
the solar mass) for C Q = 185 MeV/m" 3 , m so = 150 MeV, T c 
= 170 MeV and ji Ve = ji v = MeV. The curves a, b, c and d 
correspond to T = 0, 20, 40 and 60 MeV respectively. 

Figure 2. Same as fig. 1 but for \i Ve = MeV and {i v = 200 
MeV. 

Figure 3. Same as fig. 1 but for \i Ve = 200 MeV and {i v = 
MeV. 

Figure 4. Same as fig. 1 but for fi„ e = \i v = 200 MeV. 

Figure 5. Temperature T vs. M max /M & (where M max is the 
maximum mass) for C = 185 MeV/m~ 3 , m so = 150 MeV, T c 
= 170 MeV and ^ Ve = [i Vji = 200 MeV. 

Figure 6. Radius R vs. M/M for /i Ve = /i^ = 200 MeV and 
T = 60 MeV. Curves a, b correspond to T c = 170 MeV and two 
sets of C & m so (185, 150); (210, 100) respectively. Curve c 
corresponds to DDQM model at C = 185 MeV/m~ 3 and m so = 
150 MeV. 

Figure 7. Mass M/M vs. frequency v (KHz) for C Q = 185 
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MeVfm- 3 , m so = 150 MeV, T c = 170 MeV and fi Ve = ^ = 200 
MeV. Curves a, b correspond to fundamental mode, n = 0, at 
T = and 60 MeV respectively. Curves c, d correspond to n = 
1 mode at T = and 60 MeV respectively. 

Figure 8. Mass M/M Q vs. frequency v (KHz) for /i Ue = n v 
= 200 MeV and T = 60 MeV. Curves a, c correspond to T c 
= 170 MeV and two sets of C Q & m so (185, 150); (210, 100) 
respectively for fundamental mode. Curves b, e corresponds to 
corresponding DDQM model at C = 185 MeV/m~ 3 and m so = 
150 MeV. 
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Figure 1: Radius R in km vs. mass M/M (where M & is the solar mass) 
for C = 185 MeV/m- 3 , m so = 150 MeV, T c = 170 MeV and fi Ue = fx Vfi = 
MeV. The curves a, b, c and d correspond to T = 0, 20, 40 and 60 MeV 
respectively. 
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Figure 2: Same as fig. 1 but for fi Ue = MeV and ^ = 200 MeV. 
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Figure 3: Same as fig. 1 but for fi Ue = 200 MeV and ^ = MeV. 
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Figure 4: Same as fig. 1 but for ji Ve = pL^ = 200 MeV. 
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Figure 5: Temperature T vs. M max /M Q (where M max is the maximum mass) 



for C = 185 MeV/m" 3 , m so = 150 MeV, T c = 170 MeV and /v = ^ 



200 MeV. 



25 




0123456789 10 

R(Km) 



Figure 6: Radius R vs. M/M for fi Ue = /i v ^ = 200 MeV and T = 60 MeV. 
Curves a, b correspond to T c = 170 MeV and two sets of C Q & m so (185, 
150); (210, 100) respectively. Curve c corresponds to DDQM model at C = 
185 MeVfm- 3 and m so = 150 MeV. 
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Figure 7: Mass M/M Q vs. frequency v (KHz) for C D = 185 MeVfm 3 , m so = 
150 MeV, T c = 170 MeV and \i Vk = ^ = 200 MeV. Curves a, b correspond 
to fundamental mode, n=0 at T = and 60 MeV respectively. Curves c, d 
correspond to n=l mode at T = and 60 MeV respectively. 
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Figure 8: Mass M/M vs. frequency v (KHz) for fi Ue = /i„ M = 200 MeV and 
T = 60 MeV. Curves a, c correspond to T c = 170 MeV and two sets of C Q & 
m so (185, 150); (210, 100) respectively for fundamental mode. Curves b, e 
corresponds to corresponding DDQM model at C = 185 MeV/m~ 3 and m so 
= 150 MeV. 
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